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ABSTRACT 


The use of a magnetic field to contain plasma at extremely high 
temperatures has come to the fore in recent years as the temperatures 
in propulsion systems have risen. This thesis makes a theoretical 
investigation of the shape of a nozzle formed by a magnetic field 
acting on a perfectly conducting plasma. A description of the system 
with paeie assumptions is given, and the system is treated analytically 
in an attempt to develop some method of deriving a closed form 
solution. A numerical attack on the problem is used to find actual 
nozzle shapes for a variety of conditions of Mach number, specific 
heat ratio, current, and distance between the current-carrying 
conductor and channel wall. 

The author wishes to express his appreciation and gratitude to 
Professor Orval H. Polk of the United States Postgraduate School for 
the assistance and encouragement given him in this study. He also 
thanks Dr. Richard M. Head of the same school for arousing his 
interest in the field of magnetohydrodynamics and Dr. Burton D. Fried 
of Space Technology Laboratory, Hawthorne, California, whose work 
provided a basis for this investigation, and whose discussion with 


the author pointed out the problem areas needing investigation. 


ii 





Section 


10. 


TABLE OF CONTENTS 
Title 

Introduction 
System Description and Basic Mathematical Discussion 
Simplifying Assumptions 
Consequences of These Assumptions 
The Common Surface and its Imposed Boundary Conditions 
Current Criticality and ie Effect on Nozzle Shape 
Discussion of a Mathematical Solution 
Numerical Solutions for Specific Cases 
Conclusions 


Bibliography 


iti 


Page 


13 
15 
21 
31 


59 





Figure 


LIST OF ILLUSTRATIONS 


Diagram of Basic System 

>) 
Variation of Nozzle Shape with Change in Current 
Variation of Nozzle Shape with Change in Mach Number 


Variation of Nozzle Shape with Change in Specific 
Heat Ratio 


Variation of Nozzle Shape with Change in Conductor to 
Wall Distance, Absolute Current Constant 


Variation of Nozzle Shape-with Change in Conductor to 
Wall Distance, Number of Critical Currents Constant 


iv 


Page 


35 


36 


37 


38 


39 





TABLE OF SYMBOLS 
Sonic velocity 
Area of nozzle in equation 14; magnetic potential 
Magnetic induction 
Speed of light 
Specific heat at constant volume 


Distance from current-carrying conductor to point of given 
magnetic flux density 


Distance from current-carrying conductor to nozzle breakaway 
point 


Displacement density 


Vector operator; = = & + - V 
Electric field intensity 
y-coordinate of nozzle points 
Function defined by equation 46 
Green's function 

Function defined by equation 49 
Magnetic field strength 

Electric current 

Critical current value 

Current density 

Nozzle length 

Abscissa of nozzle breakaway point 
Mach number 

Number of critical currents; unit normal vector 


Fluid hydrodynamic pressure 





< 


“< 
.@) 


iis ea ali YM 


€ 


TABLE OF SYMBOLS 
continued 


Fluid hydrodynamic pressure in Section 8 
Pressure exerted by magnetic field 

Free stream hydrodynamic pressure 
Pressure function defined by equation 49 
Charge density 


Perpendicular distance between current-carrying conductor and 
channel wall 


Universal gas constant 

Arc length 

Entropy 

Time 

Tempe rature 

Velocity; specific volume in equation 9 
y-coordinate of fixed portion of channel wall 
Specific heat ratio 

Dirac delta function; variational symbol 
Dielectric constant 

Scalar function defined in equation 1] 
Magnetic permeability 

Fluid density 

Fluid conductivity 

Viscous stress tensor 


Nozzle shaping factor 


vi 





1. Introduction. 

Fluid flow through some constraint such as a nozzle is a phenomenon 
which is rather commonplace in these days of rockets, jet propulsion, 
and high Mach number wind tunnels. Everyone is familiar with the use 
of a physical constraint in these applications. This is generally satis- 
factory as long as the conditions involved in the flow do not exceed the 
thermal and mechanical limitations imposed by the nature of the material 
used. 

However, when the search for improved performances in these devices 
requires fluid temperatures which preclude the use of any known material, 
two approaches may be employed. New materials may be developed which 
operate satisfactorily under these conditions. While this may solve the 
problem temporarily, there still exists the possibility that the next 
step forward may create the same problem again under conditions making 
it even more difficult of solution. The other approach is to find some 
constraint which is not affected by temperatures. This could have appli- 
cation not only in the situation as outlined above, but in any set of 
circumstances where exceedingly high temperatures are encountered. 

Thus, much interest is being shown in the containment of high 
temperature matter by magnetic fields. When the fluid employed is ina 
state of ionization, that is, exists as a plasma, these techniques are 
especially ee 1iéebile. The problem of nozzling a plasma flow by means 
of a magnetic field is discussed in this thesis. 

This thesis is primarily a theoretical discussion of the problem. 

It is recognized that the conditions established would probably not be 
those of a practical situation. However, it is felt that the discussion 
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could, by extension, be applied to other situations which would be more 
closely akin to actuality. 

The thesis attempts to develop a mathematical basis for such a 
problem, using some assumptions which greatly simplify the mathematics 
while preserving the validity of the problem. A discussion of the 
analytic solution of the equations is made, which under proper conditions 
could possibly furnish at least a starting point for getting answers. 
Finally, a numerical solution of the problem, using arbitrarily chosen 
parameters is discussed. The results of these parameters are tabulated 
in detail and graphed in somewhat less detail to show qualitatively the 
effect on the system of parametric variations. These are then discussed 
and some general conclusions drawn. CGS units are used in all equations. 
2. System Description and Basic Mathematical Discussion. 

Consider a neutral, ionized fluid confined in a two-dimensional 
channel of width 2y,, as Rican in Fig. 1. An infinitely long current- 
carrying conductor is located at a distance, r, above the upper wall of 
the channel, and oriented perpendicular to the direction of flow. It 
carries a steady current, I, in one direction. A similar conductor with 
current in the opposite direction is located identically beneath the 
lower channel wall. 

The origin of a rectangular coordinate system is chosen at the 
intersection of the midplane of the channel and the line joining the 
two conductors, with x taken as positive to the right and y positive 
upward. Thus the channel mid-plane becomes the xz-coordinate plane. 
Because of the symmetry of the system about the x-axis, the shape of 
the upper half of the nozzle is the mirror image of the lower. Thus, 
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a solution for one will serve for both. The upper half is considered 
here. 


1,2, consider 


In order to mathematically describe such a system 
separately the phenomena affecting it. The electromagnetic aspect 
is at once apparent. The field effects from this source are best 


described by Maxwell's equations. These start with Ampere's law 


relating magnetic induction to the electric current producing it: 


De 
We Se Spat (1) 


Second is Faraday's Law of Induction: > 


jes / 2B 
VxE = -e Se - (2) 
The last of these equations concern the field-source relationships: 
—_> 
V-Be= O, (3) 
> 
VE = 477@. (4) 
For the electric properties of the moving medium, Ohm's law is 


needed: = B 
: > UX } 


J a EF (F t ¢ 





(5) 

In considering the flowing fluid of the system, the normal conservation 
equations of hydrodynamics must still apply, modified now to include the 
electromagnetic action upon the flow. The equation for conservation of 


mass is unaffected: 


op + V-(p¥) =O, (6) 


eller ene Cosmical Electrodynamics, Clarendon Press, Oxford, 1950 


aueyer, On Reducing Aerodynamic Heat Transfer Rates by 
Magnetohydrodynamic Techniques, IAS Preprint No. 816, 1958 


ww 


aclaieer. Magnetohydrodynamics, STL Rpt PRL-8-04, 1958 


4 





The Lorentz force and the electric field force must be included in the 
momentum arial, - 
Pe tUp = gh jxB svt (7) 
Finally, the energy equation must reflect the electrical energy which 
has invaded the fluid system, both in the form of Ohmic heating and in 
the energy ea ae "eo the stage and flow fields: 
Pp be PB -(gF+ jx 8). F4V-(F2). — @) 
Three a which are implied in the writing of these 


equations in this fashion are: 


~ 1 O"s 
% Samy 


ae 

ae 
ei * 

3. Simplifying Assumptions. 

In order to reduce this rather formidable set of equations to one 
which might be more easily handled, several simplifying assumptions are 
made. Some justification is listed for each assumption, so that it may 
be seen that solution validity is not degraded into absurdity. 

With regard to the fluid, viscous forces are considered to be many 
times smaller than the electromagnetic, and they are ignored. This 
assumption will not produce appreciable error, considering that the 
magnetic field must be strong enough to produce pressures in excess of 
the ambient pressure of the free stream. Further, the flow is considered 
isentropic with no shock effects. Since the flow is in all cases super- 


sonic, again this does not deviate greatly from the actual situation. 


jordan Electromagnetic Waves and Radiating Systems, Prentice-Hall, 
New York, 1950 





No flow exists in the z-direction, the direction of current flow. 
Since the main interest in this thesis is in the field-flow interre- 
lationships, rather than in the flow itself, the one-dimensional channel 
flow approximation is used to describe the two-dimensional flow. This 
is a standard fluid mechanics peenniare: of sufficient accuracy to 
justify its use in most nozzle problems. A steady state has been 
reached in the flow. Finally, the plasma is assumed to behave as a 
perfect gas, with state and entropy equations SE. 
pv = RT (9) 

S: Cy In En (10) 

Electrically, the fluid is assumed to have infinite conductivity, 
since the conductivity of a plasma is extremely high. 
4. Consequences of These Assumptions. 

The assumption of infinite conductivity within the plasma tre- 
mendously simplifies the system of equations 1 through 8 since the 
currents induced in the conductor through the action of the external 
field will just balance this field, thus excluding the magnetic vector B 
from the inside of the fluid. This is a restatement of the fact that 
a perfect conductor is likewise a perfect diamagnet. It is quickly 
arown’ from a consideration of equations 1 through 8 at steady state, 


eT itoman and Puckett, Aerodynamics of a Compressible Fluid, Wiley 
and Sons, New York, 1947 


Bk iStee Kinney and Stuart, Principles of Engineering Thermodynamics, 
Wiley and Sons, New York, 1954 : 
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with @ =00 and q = 0, conditions implied by the property of infinite 
conductivity. Then equation 5 becomes 
= , 
FUK2-0 
— > —— 


or VX oe ie. 
> 
E 


But, within the conductor, does not exist, and therefore 


> 


rxB=O. 

This implies that these vectors are parallel or that 

_ 

B= Apr, (11) 
where A is a scalar point function defined by this equation. 

It is further known from equation 3 that 
—_ 
V-B=0, (3) 
and that at steady state, equation 6 becomes 
V- (pr: 0. (12) 
Substituting 11 into 3 and expanding 
U-Apr AV: PR) + pu-vA. 
But, from 12 V- (pF): 0” 
and thus Ved Fs pv-Vh. 
From 3 ie = Vid py =O. 
> 
Therefore, ae -VA=O ’ 
indicating that along any streamline, » is a constant. But, since 
e@epprogches zero as x becomes infinite, mn must be zero on all stream- 
lines at infinity. » is constant along each streamline and is thus 
zero everywhere in the fluid. From Li, ae therefore be zero 
throughout the fluid. 
The use of the steady-state, isentropic, non-viscous, one-dimensional 

channel flow approximation to describe the fluid flow reduces the 
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hydrodynamic equations to the simplified aa Ss 


(mass ) i A Ve = fe Ao K, (14) 


(momentum) pe gy + Vp = yx c (15) 
(energy) xu*+ 2 = pals +t Be (16) 


The assumptions used have thus reduced the problem to a consider- 
ation of two coupled free-boundary problems with the only common 
element provided by the surface between them. These problems are 
essentially the existence of a magnetic field due to a current I, in 
the presence of a perfect conductor, and the isentropic flow of a non- 
viscous fluid through a constraining channel. 

The common surface is an unknown and will be represented as 

y =f (x) 
or simply, for brevity 

Veron (17) 
5. The Common Surfage and its Imposed Boundary Conditions. 

The surface, y = £, can be considered in two sections. When 
{x{€ m as shown in Fig. 1, y = f£ will actually be the locus of all 

points at which the fluid pressure due to hydrodynamic flow exactly 
balances the pressure externally applied to the plasma by the magnetic 
field. When [|[|x|>m, y = ee 

This pressure exerted by the magnetic field may be derived from 


; ; ya ; 9 
a consideration of the simplified equations at zero velocity 


Liepman and Roshko, Elements of Gas Dynamics, Wiley and Sons, 
New York, 1957 


aon Particle Transport, Electric Currents, and Pressure Balance 
in a Magnetically Immobilized Plasma, Phys. Rev. (97), 1443, Jan-Mar 1955 
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Equation 15 becomes 
es. 
=feko os (18) 
Vp =)*-€ 
Under the steady state assumption, with a direct current, equation 1 
is now 


Vi a (19) 


Substituting from 19 for j in 18, 5 
| é 
NV P=a7 (VxB)x 
—> 
¥ 471 Vp = (VxB)x 


Qov 4 |Goy 


(20) 


> > 7. 2 
Expanding, ia Vp = (B-V)B ae VB (21) 
But the gradient of the field is small in the direction of the field 
and the first term may be neglected. pener one, 
Vp= ~ er VB" 
and Jo= Or (22) 
This states the pressure exerted by the magnetic field on the fluid 
in terms of the field intensity. 
The relationship of B to the current producing the field is given 
in any basic text on Bilecucoma snetice. and in CGS units is 
B = ut | (23 
The hydrodynamic pressure at any point in the flow must then be 
determined. Starting with the energy relationship, equation 16, an 
equation can be developed which expresses the pressure ratio at any 


point in a nozzle as an implicit function of the ratio of the nozzle 


area at that point to the free stream area. If a unit depth in the 


ee Gendan: Op. eLt. 





z-direction is considered, this area ratio will reduce to a ratio of 
nozzle y-coordinates. 


If equation 16, 


y _/ 
zu veg, Boy mB, . 


is divided through by the final term, the result is 


vet RP Rt! ~ 


If the sonic velocity ieee a 


» 8 

af. Be, (25) 

is substitued in equation 24, it becomes 
Ory? b B_ delve 
2G * php -~ 2 a2 


If the first term is multiplied by Vo /i, and the definition of the 


+/, (26) 


Mach number, 
Vo 


is substituted, the equation becomes 


G-/ y* M +P = oy Hf (28) 
Z Vz» 
For a unit depth in the z- ae the mass equation, 14, may be 


— = eee 09 
Vo PF ey 


with £ defined in equation 17. Substituting this, and making use of the 


written as 


isentropic relationship, 


-Y 
-(£)° (20) 


to remove the density term, a final pel is 


(4058) LB) ee, (31) 


“oot. & and Roshko, op. cit. 
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It is noted that, if incompressible flow should be under consider- 
ation, the derivation is much simpler. In this case, the energy 
equation reduces to Bernoulli's equation, 

yi. 2 _ / 2 
przfre=pots Pre . — 
Dividing through by Po» and again using equation 14 to introduce the 
nozzle shape and eliminate the velocity, 
e /+ ZX, (y2-v) 
& 
2 ce 
} ft £v2 (/- 2) (33) 
fe *' Vs 
2 
| ha ay AP L1- (9) 
2 fo 


This is not only a simpler derivation, but because of the explicit 


t4 


relationship between pressure and nozzle shape, would probably simplify 
some of the later discussion. 

The boundary conditions imposed by this problem must now be 
considered. Since the fluid is confined everywhere by the surface of 
the nozzle, 

> 

i Ae O (34) 
everywhere on y = £. Further, since there is no magnetic field within 
the fluid, 

> > 

B-n=0O (35) 
along the same surface. Finally, since the magnetic field is shaping 
the nozzle in the region where -m ¢x m (See Fig. 1), the magnetic 


pressure must equal the fluid pressure at every point in this region in 


order to establish an equilibrium surface. Thus, 


2. 
i. i (36) 
817 p 

at every point on the surface where -m < xx m. 


11 





In the discussion of the boundary conditions and in the further 
mathematical development, it is convenient to describe the magnetic 
—_> 
field by means of a vector magnetic potential A ne such that 
—_ — 
B Seay XA. (37) 
—> 
The symmetry of the problem indicates that A will have only one non- 
a 
zero component, A? permitting it to be treated as a scalar quantity. 
This potential must satisfy Laplace's equation outside the source wire 
and Poisson's equation within it. Thus, outside, 
VY A=0 (38) 
and inside, 
2 4/T 
VA=- = Lf (39) 
Introducing a Dirac delta function, defined by 
d(4) =O, a 
J Sladz=/ 
- 00 
these two expressions may be written as a single term 
2 diy 
=e - f 40 
VA=- TZ Sle) SL y- ly tr (40) 
Relating the boundary conditions expressed in equations 35 and 36 
to this magnetic potential, it is seen that they are satisfied if 
A =0 
on y = f or, expressed in another notation, 


ALY fF] = 0; (41) 


and, if 
IVAI* = S77 p (42) 


on y = f, “m xX ¢ Mm, 


So ondane OP. Cit: 
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The boundary conditions are complicated by the fact that this is a 
free boundary problem; that is, the boundary is unknown, and in fact, it 
is the objective of this thesis to determine this boundary for a given 
set of conditions. 

There are, however, some general restrictions which may be placed 
on the boundary which will aid in visualizing a solution. The symmetry 
of the nozzle about the x-axis has already been mentioned. It is further 
apparent from Fig. 1 that, since the process is assumed isentropic, the 
nozzle is likewise symmetrical about the y-axis. The nozzle cannot have 
a sharp corner at x = +m, since this would imply the existence of an 
infinite magnetic pressure at that point, a point patently inconsistent 
with the equality of magnetic and hydrodynamic pressures stated in 
equation 36. Thus the slope of the nozzle surface must be zero at 
x =im. This slope must likewise, from symmetry considerations, be 
zero at x =O. It is then apparent from the mean value theorem of the 
calculus that inflection points exist between x = O and x = £m, and 
that the slope of the nozzle is a maximum at these points. 

6. Current Criticality and its Effect on Nozzle Shape. 

From a further consideration of the implications of the equality of 
magnetic and hydrodynamic pressures at all points along the nozzle, it 
can be seen that some minimum critical value of current must be used or 
the magnetic pressure even at the closest point of flow to current will 
not be sufficient to permit this equality to be established. For any 
current I through the wire, the maximum field intensity along the channel 
wall will occur at the point (0, y J» and can be expressed from equation 23 

41 


as B = ra - (43) 
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If I is taken as zero, the plasma flow would proceed as though it were 
ordinary hydrodynamic flow through a constant area channel. As I 
increases, B will likewise increase, but the flow will remain unaffected 
until the magnetic pressure becomes equal to Po: Thus I must have at 


least a value so that 


loLe 


Bi7ctr> _, 
lizcr oe, (44) 


: Cc 


This is referred to as the critical current. For values of I greater 
than I. » the flow will be nozzled down as much as is required to bring 
the magnetic and hydrodynamic pressures into equilibrium. When equili- 
brium is reached the nozzle shape is established. 

It is further apparent that, for any current value in excess of 
the critical, the magnetic and free stream hydrodynamic pressure will 
be equal at the breakaway points, x tm . Representing the distance 


from the wire to the rat by by d ae this may be stated as 


bd - = Po, 
from which 
ee EZ 2 Tia 
m ~ I¥c* Pe ° 


From the Pythagorean theorem, it is also true that 
2 
= 2 
lo ae a tin* 
Thus, 2 
ay 
Tc*Pe 


But any current in excess of the critical value may be represented as 


-~pé 





m* = 


some multiple, n, timesthe critical where n > 1. Thus, 


pp Zee z 


77<* Po ee = 


14 





Substituting from equation 44, 
yn*= Lift Ups 
Z7c*Po 


MW.) m= (45) 


get 


Hence, the nozzle length is directly connected to the relationship 
between current used and the critical current for the system. 
7. Discussion of a Mathematical Solution. 

The task of solving these equations is a formidable, if not 
impossible, one for the general case. Hence, an analytic solution in 
closed form has not been obtained, but an approach which may be useful 
in finding the solution in certain specific cases is given. 

This involves the variational principle of the calculus of 
variations. Dr. Fried uses ee in a similar case where he has employed 
a non-dimensional magnetic potential. 

Let the right hand side of equation 40 be represented by a single 
symbol, say g. 

g= - HT Seas Sly-( grr): (46) 
Then A = g. (47) 
An energy integral for the system may be defined for any pair of 


functions, A(x, y) and f(x), as, 


aes 2 gc a 
E (Af) = fakdy (s-IVAItE Atl oe Al] wa 
/ K YY = Co 
where R is the region of the xy-plane above the curve y = f£, and where 
h \é (x) | is a function of the pressure at any point on the nozzle 
satisfying the =e as 
(49) 


cre ~ Fe b, 


p being developed from equation 31. This term will be defined as P/2. 


eriedl op. cit. 
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2 
The first term in equation 48 is an integral of a7 over the 


region where it is not zero, while the latter term integrates p through- 
out the volume occupied by the fluid, thus comprising an energy integral 
for the system. 
From equation 41, it is known that A and f, although both are 
arbitrary, must satisfy 
Aly f]-0o. (41) 
Then small variations, 5A and $f » give a ct faape. relationship, 
i aa cg be Aly £2 A) 
SEs Joly dy (sp VA-V6A+ £2 A) SG + 
Oo 
-/ du bFF . (50) 
~ ©o 


These terms represent respectively the variation of the first term 
of the energy integral with variation in A; the variation of the same 
term with variation in f£; and the variation of the last term of the 
energy integral with £ only, since it is independent of A. The negative 
sign of the last two terms derives from a decrease in the value of the 
integral concerned when a positive (increasing) variation is permitted 
in £. The second integral in equation 50, although coming from the 
area term of the original energy integral, is now integrated with respect 
to x only, since y, which is equal to f, has now been replaced in the 
increment by Kee 

If Sf is considered to mean small variations in the surface bound- 
ing the fluid flow, it will vanish when f coincides with ye in the fixed 
portion of the system. As a result, the last two integrals need be 
evaluated only from -m to m, since they will be zero outside this 


range. 
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The first integral in equation 50 may be further treated by 


applying the vector identity'**!> 


VA-WSA = V-(SAVA)- SAYA cn 
From this, 


Jd dy (3372 VA-V5A + $2 iA): fees GS [¥. (6AVA) 
SAYA] +3 + SAP (52) 
Combining terms, expanding, treating a as a constant value, and 


applying the divergence theorem to the result, it is further reduced to 


2 vy oe 
Lay ay [fA $2 - £ VA)S- inl ds (7. VASA (53) 
where ds is the element of arc length along y = f, and n is the unit 
normal to s inward on R. 
Again recalling that, from equation 41, 
aly, f] = 0 . (41) 
it is seen that along y = f, 
Sar oe ° Sf 0, 
or that 
-dA= de r ot. (54) 
Considering the second integral a 7 53, it is noted that the 
term in parenthesis may be rewritten as 


n-VA= oo 


Substituting these two into the second integral, and recalling the 


(55) 


remarks on the limite of integration above, it may now be written 
m 


eeidaee Advanced Engineering Mathematics, McGraw-Hill, New York, 1951 


+ BeaTid Vector and Tensor Analysis, Wiley and Sons, New York, 1947 
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With the expressions developed in equations 53 and 56 substituted 
into equation 50, the variational form of the energy integral, it can 


now be rewritten as 


§F: a ernie ae [3h i 53 F 


CK 
— £,/VAl*- so A a £ | (57) 
Since A[x, f | = 0, the _ = A must be in the direction of 
the normal to y = £, and therefore | 
n. WA = 2 =|VAl. (58) 
Considering 
dA ds _ 9A ds 
af dk dy ok 
_ OA Myr +y)~ 
dy Att 
- _| (gA¥{@«era@y) 
gay] "| 
__[ (dae, (dA\* 
7 oe) ul te 
- [VAI a 


) 


it is seen that 


S ds = [VAIS (60) 
g:- OT §() 5(y Ger) (46) 


vanishes on y = f£, the end result is 2 
SE+ fav dy a(S - 5, 0'A) 4/"ebe foo | VAL” 
: Pt. (61) 


This will vanish, satisfying the variational principle, if and only if 


SIS 


Now, since 


A and £ satisfy the equations 41 and 42. 
These results may be further applied to this problem by noting from 
equation 41 that A is a functional of £; that is, A depends ina unique 
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manner on f. Then, if f were chosen on any convenient basis, A would 
be determined by equations 41 and 61. 

Consequently, the energy integral becomes a functional of f alone, 
rather than of f and A, and setting 

sF& =O (62) 

as prescribed in the variational principle, should develop an equation 
which is in terms of f, and which would be satisfied by the f£ which 
together with its associated A, satisfies equations 41 and 42. 

At this point, consider the introduction of a Green's function of 


two arbitrary position vectors of the region R, 


G(r, rAEFs, In (7-7) . (63) 
VG(r, n)= &(4,-fA), (64) 


it is possible to use a function 


g(7)=- 216 Ff) 65) 


Because 


where 
a 
2) eee) (66) 
as a solution for 
WA(r) = g(7)=- SISGF-Z), (67) 
which restates in vector form equation 40, since the Dirac delta function 
will be non-zero only at the position of the current carrying conductor. 
Reconsider now the statement of equation 48. Since the second 
term of the integral is already a function of f alone, only the first 
term will be discussed at this point. 
This term may now be written, considering the area summation as 
the summing of all position vectors, fF, » as 


Lar 55/421 A), (68) 
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However, employing the vector identity 


V-(AVA = |VAI+AVA, (69) 
this may be expandee to 
fai ($l (AVA) - AVA] + OF g A). (70) 


But, from equation 46, 
wv A = 8 .- (46) 
Collecting terms, and noting that, from the symmetry of the problem, 
the integral of the first term vanishes, the integral becomes 
iY ff oe = == 5) / 
Lo, ( RA HG A) 4t (/ &I io; 9A (71) 
But because of the definition of D (r) above, 
2 
9= WV", (72) 
Na; - & = 
18 (/ Z) = k (73) 


and setting 


the integral now becomes 
KJ de Avg. on 
If the expression 
SVAtVA-V¢ 
be added and subtracted from the integrand, the integral develops into 
k fan (9v'A tVA-VE AVP - VA VP- BVA) 025) 
This may be manipulated into 
kfdr[9v7Aar VAVG -T(GVAY 
= k foi [PVA +¥-(AVO- p val (76) 
The integration of the first term gives 
-k7 $&). (77) 
The second term vanishes over the area integration, and the third term 
may be transformed by means of the divergence eimorén to 


k [ds g gh (78) 
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by use of equation 58 and recalling that n has been defined as a unit 
normal inward on R. 

If ds is expressed in terms of dx, the first term of equation 48 

becomes o ae 
E=-k gai rf dvfk Li-(ZYT" o 34 (79) 
Combining this with the second term of equation 48, the energy integral 

may be written no fe I, 
E= -kt gli) + Leale {kt 7 a) J ? 34 +hf (80) 

Although the first term is infinite in the limit, it is independent 
of £ and thus, when the variational principle is applied, may be ignored 
as it will contribute nothing to Se. Thus, formally, the problem has 
been reduced to a dependence on f alone, and when $E is set equal to 
zero, an exact equation for f will emerge. Practically, however, the 
dependence of A on f, as defined by equation 61, is so complex in 
general that the expression for SE cannot be written. 

8. Numerical Solutions for Specific Cases. 

Because of the complexity of the relationships necessary for 
analytic solution, the problem was attacked numerically, using the 
Model NCR-102A computer facility available at the U. S. Naval Postgraduate 
School . 

In this attack, equation 31 was first solved. There are four 
variables in this equation. The specific heat ratio and Mach number 
were arbitrarily fixed at one of five preselected values. Mach numbers 
began at 5.0 and proceeded in steps of five to a value of 25.0. Values 
of the specific heat ratio ran from 1.30 to 1.50 in increments of 0.05. 
All possible combinations of these two parameters were eventually used 


Pad) 





in the solution. This reduced the equation to one in which only two 
variables exist for any specific case. Because of its exponents, the 
pressure ratio was selected as the independent variable, with values 

from 1.000 to 5.000. This permitted solution for the y-coordinate ratio 
corresponding to a given pressure ratio. An interpolation routine in 
the computer then converted the ratio of y-coordinates to the independent 
variable, tabling the pressure ratios corresponding to preselected 

values of this parameter. This approach was used because the use of 
pressure ratio as the independent variable permitted a simpler solu- 

tion of equation 31. The final results were given for preselected 

values of the y-coordinate of the nozzle, allowing simpler graphing, 

and, to the author, permitting a more readily understood and quickly 
prepared tabulation of nozzle shape data. The values of pressure ratio 
corresponding to the nozzle y-coordinate for the combinations of M and 

% are given in Table I. 

With this accomplished, for a combination of Mach number and 
specific heat ratio chosen from those used, the fluid pressure at any 
point along the nozzle can be determined in terms of the free stream 
pressure if the y-coordinate of the nozzle at that point as compared 
to Y, is known. Since, at every point along the nozzle, the magnetic 
pressure must equal the fluid pressure, the magnetic pressure is known 
in its relationship to the free stream pressure once a y-coordinate 
is picked. 


The magnetic pressure is related to the current in the conductor 


through equations 22 and 23, 


= 
Pm = pg (22) 


= (23) 
0 cd 
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in the solution. This reduced the equation to one in which only two 
variables exist for any specific case. Because of its exponents, the 
pressure ratio was selected as the independent variable, with values 

from 1.000 to 5.000. This permitted solution for the y-coordinate ratio 
corresponding to a given pressure ratio. An interpolation routine in 
the computer then converted the ratio of y-coordinates to the independent 
variable, tabling the pressure ratios corresponding to preselected 

values of this parameter. This approach was used because the use of 
pressure ratio as the independent variable permitted a simpler solu- 

tion of equation 31. The final results were given for preselected 
values of the y-coordinate of the nozzle, allowing simpler graphing, 

and, to the author, permitting a more readily understood and quickly 
prepared tabulation of nozzle shape data. The values of pressure ratio 
corresponding to the nozzle y-coordinate for the combinations of M and 

% are given in Table 1. 

With this accomplished, for a combination of Mach number and 
specific heat ratio chosen from those used, the fluid pressure at any 
point along the nozzle can be determined in terms of the free stream 
pressure if the y-coordinate of the nozzle at that point as compared 
to a is known. Since, at every point along the nozzle, the magnetic 
pressure must equal the fluid pressure, the magnetic pressure is known 
in its relationship to the free stream pressure once a y-coordinate 
is picked. 

The magnetic pressure is related to the current in the conductor 
through equations 22 and 23, 

Nee 
bm = == (22) 
Ba a. (23) 





Therefore, if a known current flows in the conductor, the distance from 
the conductor at which the magnetic pressure is equal to the stream 
pressure for a certain nozzle y-coordinate may be calculated. If the 
conductor is located a specified distance from the outer wall of the 
channel, then Bl) conditions are known to completely determine the 
nozzle shape. 

Thus, to develop the resulting nozzle shapes for specific conditions, 
current flowing in the conductor must be known. As described above, a 
certain minimum current must flow, else the fluid is undisturbed. This 
critical current was shown in equation 44 to be 

Lo=cr yf &. (44) 
If the current is chosen in multiples of this critical current, as was 
done for this thesis, three advantages accrue. First, all troublesome 
constants disappear. Secondly, the equations reduce to simple expressions 
which are ideally suited for digital computer work. Finally, the work 
is done non-dimensionally and thus wider applicability of the results 
are gained. The currents used here are respectively, 21 31 4I oy 
and oI. 

The final item which affects a specific nozzle shape is the distance 
r, the separation between the current carrying conductor and the channel 
wall. This influences the magnetic pressure, since the absolute value 
of the critical current depends on this distance. It also has a more 
direct influence on the nozzle shape since the abscissa of any point on 
the nozzle will be the third side of a right triangle whose hypoteneuse, 
d, and side, r, are known. The advantages of non-dimensionality are 
again applicable and, therefore, in this problem, the r's used are 
¥/5: y/4) and y ,/3- 
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A variation of Pp, was not considered in this thesis. Since the 
currents were used in multiples of the critical cunteort which is a 
function of Pe the variation of current parallels a variation of 
pressure. Thus, if current were held constant and pressure changed by 
a factor, the number of critical currents represented would be changed 
by the inverse square root of that factor, and the new nozzle would be 
shaped exactly as a nozzle whose forming current was represented by 
this number of critical currents with the initial pressure. 


As a corollary to this, the shaping factor of the nozzle may be 


expressed as the ratio, 





VE 
Te 


If this factor is held constant, the nozzle shape for a given Mach 


W = (81) 
number, specific heat ratio, and conductor-wall separation will not 
change. 

By studying the effects of varying each parameter in turn for 
constant combinations of the others, both general qualitative statements 
of effect on nozzle shape and the actual shapes for the given conditions 
can be realized. 

In arranging the work for digital computation, some multiple of 


the critical current was assumed to flow. 


(ewe ger 7 (82) 


Therefore, 
ffs 
B= éq/ (23) 
and 





2 
fm = 3B a (22) 
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le of * 
2 ae 

BE = Gy oo 


But from the boundary conditions, at any point along the nozzle, 


Pm = Pf (84) 
Ln . ff 
Po Po 


or 


or 


Therefore, 7 
Ay = Be os 
ga... (86) 
This equation is the essence of the eo ER ak calculation of a 
nozzle shape for a particular Mach number, specific heat ratio, current, 
and distance of conductor from the channel wall. 
For given values of M and 2 , let a value of (f) be selected. 
This accomplishes two things. The value, (FE) , is fixed, and a straight 
© 
line locus of points in the fluid at which that pressure ratio may exist 
is established. With a fixed pressure ratio and current, the distance 
at which a magnetic pressure equal to the fluid pressure may exist can 
be calculated in terms of r. Since r is known, a circle of now known 
radius is the locus of points at which the equal magnetic pressure 
exists. The intersections of these two loci establish the points 
on the nozzle. 
To calculate the nozzles for specific conditions, the values of 
( £4) previously tabulated for those conditions were fed to the computer, 
along with the number of critical currents in each of the currents 
used. The indicated division was performed and the (a) ratio was 
obtained. These ratios are tabulated in Tables II through VI for each 
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of the conditions outlined. With these values, graphs of the nozzle 
shapes were made by constructing the loci as described above. Graphs 
were not made for each of the sets tabulated. This was unnecessary 
since those graphs shown indicate the effect of a change in one parameter 
with the others constant. This effect is similar, as shown by the values 
tabulated, no matter what values are assigned to the fixed parameters. 
Each nozzle has a minimum diameter. This occurs when the circular 
locus of equal magnetic pressure is just tangent to the straight line. 
That such a point exists is evident from the following consideration. 
As the y-coordinate of the nozzle decreases, the fluid pressure increases. 
Consequently, the magnetic pressure must likewise increase to maintain 
equilibrium. This demands a smaller diameter for the circular locus. 
Thus a point will be reached where the diameter of the circular locus 
is just equal to the sum of the conductor-wall separation and the fluid 
boundary displacement from the wall, and here the nozzle minimum will 
occur. The (d/r) value of these minima are included in Tables II 
through VI, for values of r = 0.2y, for all conditions, and r = 0.25y 
and 0.33y for all ‘S a 1.30. 
In actually determining the point at which the minimum occurred, 
it was helpful to construct a table giving the vertical distance from 
the conductor to each of the y-coordinates used in the problem. If the 
value of (2) as calculated above exceeded this distance for a specific 
y-coordinate, then the two loci intersected. If the reverse were true, 
there was no intersection, and hence the nozzle surface did not exist at 
that y-coordinate. A comparison of this table with the tables of (2) 
constructed for the various conditions gave the crossover values, and 
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a cross interpolation provided the actual minimum point. This table 
is Table VII. 

It may also be noted from equation 86 that at the breakaway point 
in each case, the ratio, ee » is exactly equal to the number of 
critical currents, n, since at that point, the ratio, (EE), 4s unity. 

. Consider now the qualitative effects of variation of these para- 
meters singly with all others held constant. 

As seen from Fig. 2, or by noting the values in Tables II through 
VI, an increase in current invariably deepens and lengthens the nozzle. 
This is predictable, since the magnetic pressure is directly proportional 
to the square of the current, and thus, at any point, the fluid must be 
nozzled down to ieremre the hydrodynamic pressure for equilibrium. 

As seen from equation 86, the distance from the conductor to the 


v 


breakaway point is always 
d,=nr 
where n is the number of critical currents flowing in the conductor, 
n>| . Equation 45 also points out that the abscissa of either break- 
away point is, in magnitude 
(n»-i)" r, 
If the length of the nozzle is defined as the distance between breakaway 
points, the length in this case is 
l, = 7p saree (87) 
If the current is changed by a factor & , the number of critical currents 


is likewise changed by &, and the new nozzle length would therefore be 


| > 2r 1an-/. (88) 
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Thus for a change of current by a factor & , the nozzle length is changed 


by a factor, 
2p2 


— = ——$_——_ 


Qa-T 
indicating the change in nowt length depends not only on the change 
in current, but also on the amount of current initially flowing. 

Defining nozzle depth as the distance, yo - f£, at the throat, the 
change of depth lends itself to no such simple prediction formula 
because of the complex relationship between the nozzle y-coordinate 
ratio and the pressure ratio expressed in equation 31, and is best left 
to a numerical calculation in each case. 

If Mach number of the flowing fluid is the variable, some interest- 
ing results occur. First, the breakaway point is located at the same 
spot, regardless of Mach number, since this parameter can influence the 
radius of the circular locus only through the pressure ratio, and this 
is unity at the breakaway point, regardless of Mach number. As Mach 
number increases, the nozzle depth likewise increases, because with an 
increase in Mach number, the ratio of y-coordinates necessary to achieve 
a given pressure ratio decreases, as may be seen by a comparison of the 
data in Table I for various Mach numbers. Again, the complexity of 
equation 31 precludes setting up a simple expression for predicting 
depth changes. The most interesting aspect of variation of Mach number, 
however, lies in the fact that changes in Mach number for the lower 
values used cause only small disturbances in the nozzle shape and changes 
occurring above the Mach of 15 cause almost no variation, as shown in 


Table I. Fig. 3 shows the result of this variation. 
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When Ee, the specific heat ratio of the flowing fluid, is changed, 
the breakaway point again remains unaffected for the same reasons 
discussed in the case of Mach number variation. As seen in Table I by 
comparing the pressure ratios for a given y-coordinate ratio at various 
values of specific heat ratio, an increase in ® necessitates a larger 
y-coprdinate to realize a given pressure ratio. As a result, the depth 
of the nozzle decreases with an increase in & . This change in depth 
is small but noticeable. Again it is unpredictable. The effect of 
changes in , are shown in Fig. 4. 

If the distance r between the conductor and the channel wall is 
changed, with the absolute magnitude of current constant, a double 
effect is introduced. Since the value of the critical current is 
dependent on r, the number of critical currents flowing varies inversely 
with the change in r. And since r has been changed, the value of d, 
must change, since the (4) ratio will still have the same value. 


Thus, if r is varied by a factor B, 
l= Bro 


L =f Cr, > Lhe (90) 


AT 
Applying equation 86, 


qd a nlp 
7 (91) 
i PR, 
and at the breakaway point, 


-2;,:n7 (92) 
d= Ly, oO 





Thus, an increase in r causes a decrease in nozzle length since the 
distance from the conductor to the breakaway point remains constant, 
but the side r of the right triangle has increased, necessitating a 
decrease in m. 
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The length of the nozzle is now 


Lelia = 2 ae 
= 2r1n-f* . (93) 


Since 
(= 2r%n® 1 ; (94) 

it is quickly seen that a change in r by a factor G causes a change in 
nozzle length by a factor of 

f = hep ; | (93) 
This is in accord with the conclusion mentioned earlier, since if fe 
is greater than unity the numerator of this expression is less than the 
denominator. This demonstrates the inverse relationship between 
conductor-wall separation and nozzle length. It is also to be noted that 
the value of 6 may not exceed that of n. When these two are just equal, 
the current flowing in the conductor will be exactly one critical 
current. Any increase in p beyond this value will result in an 
undisturbed flow. 

The nozzle depth is likewise diminished when r is increased. 
Consider that the y-coordinate of the nozzle throat remains constant 
when r is increased. The magnetic pressure at that point will drop, 
since the distance from the conductor to that point has increased. But, 
for fixed Mach number and specific heat ratio, the hydrodynamic pressure 
is a function of the y-coordinate of the nozzle only. Therefore, to 
establish equilibrium, the throat must expand to reduce the hydrodynamic 
pressure, and the depth of the nozzle is diminished. These results are 


illustrated in Fig. 5. 
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If, in lieu of holding the absolute value of the current constant, 
the number of critical currents is made the constant factor, just the 
reverse occurs. If r changes by a factor, @ , so that 

GEG, (94) 
the distance between the conductor and the breakaway point will still 
be, from equation 86, 

a 


=> = 86 
2 =n, (86) 


since n is a constant. The nozzle length is then 

l = Lr,Atn*-/. (95) 
Since the initial length was 

l= 24! (96) 


it is quickly seen that 


[> 


G ss U -_ 

oS re Q, (97) 
7] 

Thus the nozzle length is increased by the same factor as the distance r. 

Nozzle depth likewise increases. This can be seen from a consider- 


ation of the fact that the absolute current has increased by the same 


factor as the distance r, while the distance d, which is the sum of(y -£) 


and r, is increased by a smaller factor, since ye is a constant. Thus, 
since the magnetic pressure is proportional directly to the saunte of 
the current, and inversely to the square of the distance d, it is seen 
that the magnetic pressure at the position of the original throat has 
increased. The nozzle throat must then contract to increase the hydro- 
dynamic pressure and establish equilibrium again. Fig. 6 is a graphical 
presentation of this effect. 
9. Conclusions. 

As a result of the solutions assembled and dis cudseduin Section 8, 
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it is concluded that, in a magnetic nozzle of the type described 

herein, the most effective method of changing nozzle shape is by varying 
the current. Variation in Mach number of the flowing plasma or its 
specific heat ratio has too little effect to be of much use. Nozzle 
length is completely unaffected. Beyond a Mach number of 15, the nozzle 
shape becomes nearly independent of the value of this parameter. These 
would also be considerably more difficult to control than the current. 
The variation of distance would probably be limited by the physical 

Size of the system. While it is actually the ratio (z) which is the 
nozzle shaping factor, P, would be much more Cae vary effectively 
than I. For these reasons, current variation is the most effective 
single factor in shaping the nozzle. 

It is further concluded that for extensive work in this field, much 
more comprehensive tables will be needed. While the values taken herein 
provide indications of trends in the system, they would hardly suffice 
for much work where the parameters differed from those used in this 
thesis. Because of all the non-linearities involved, linear interpola- 
tion between tables would probably not be satisfactory, and extrapola- 
tion would be meaningless. For each new combination of parameters 
involved, the complete solution would have to be made, as no solution 
in closed form was developed. 

Mention is made of two mathematical methods which might prove useful 
in further attacks on this type, of problem. While both of these were 
considered:-by the author in his researches, lack of time prevented 


extensive development of them. The first is the method of characteristics 
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which is expanded by touennendees based on the work done by Courant, 
Issacson and Rea The numerical work in this scheme is very involved, 
and a discussion of it with the director of the computer facility used 
convinced the author that the computations would be impossible to do in 
the time available. 

The second is the method of the hypercircle, developed by warn. 
which has proved a powerful tool in the solution of types of differential 
equations which were insoluble by other means. Again, the complexity of 
the calculations deterred the author from pursuing this path. 

Finally, it must be mentioned that in the same case with plasma 
flow velocity zero,’ an analytic solution for the shape of the nozzle 
has been made by poeta using the technique of conformal mapping as 
suggested by eae In this solution, the free boundary can be mapped 
into a fixed boundary, since, when the velocity is zero, the pressure 
throughout the plasma is constant. This technique does not improve 
the situation when the velocity is not zero, since the boundary after 
mapping is also still unknown. Another interesting aspect of the zero 


velocity case is the study of the magnetic piston. Under this velocity 


mvoughhead, Solution of Problems Involving the Hydromagnetic Flow 


of Compressible Ionized Fluids, Phys. Rev. (99) 1678, Jul-Sep 1955 


concent . Issacson, and Rees, On the Solution of Nonlinear Hyperbolic 
Differential Equations by Finite Differences, Comm. P. Appl. Math. (5) 
243, 1952 


acy dae The Hypercircle in Mathematical Physics, Cambridge Press, 
Cambridge, 1957 
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Sab). Hydrodynamics, University Press, Cambridge, 1930 
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condition, the throat diameter can be made zero by a magnetic field of 
sufficient intensity. 

Further studies in this field that suggest themselves include the 
extension to a flow which is bounded in the z-direction and which is 
nozzled by the current flowing in a circular coil about the fluid. This 
could be further extended to current flowing in a solenoid. The problem 
of nozzle shape when fluid conductivity is not infinite is a more 


practical, but also a considerably more difficult one. 
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TABLE II 


y poeta ee M 
ee 





50 TTC ~ 1.0945] 1.6 889 | 2.7361 — 
42 1. 1. 2. 27 1 1 Ap 2.8253 
44 i 1. 1 25 2. 1 1 2 2.9131 
46 /1. 1 rae 2. 1 l. 23 3.0001 
48 il. 1 Ze 3. 1 1, 2. 3.0853 
290 Il. 1 2. 37, 1 Le 2. 3.1680 
m2. | 1. 1 Zs 3. 1 Le 2. 3.2514 
: 4 | 1. 1 2. Di 1 Li Ze 3.3329 
| .56 il. 2 Ze 3. 1 2. 25 3.4130 
.58 11. 2 ae 3. 1 2. Zz: 3.4926 
.60 j 1. 2 2 3e 1 Die 2: CS SIA 
-62 11. 2 td 3. l 2; La 3.6515 
.64 J 1. 2 2: ie 1 fae Ze 3.7294 
.66 j1. 2 ae On l 2. oi 3.8046 
.68 11. 2 3 a5 1 Ze 3. 3.8796 
O11, 2 Si: 3% l age 3% 3.9554 
a72 1d. 2 ce oe 1 yas a7 4.0277 
74 41. 2 3. 4, 1 Z. 3% 4.1014 
76 31. 2 a 4. 1 Zs oe 4.1729 
.78 1. 2: 3. 4. 1 2% oe 4.2459 
.80 }1. 2. 3% 4, 1 Zz oe 4.3159 
.82 1. Z. 5% 4, 1 Za Sy 4.3886 
.84 jl. 2. a. 4, l Die a 4.4559 
.86 1. Zs 3. 4, 1 Zz 5 4.5265 
.88 | 1. 2. 3. 4. 1 De or 4.5952 
.90 j1. 2. 3. 4. 1 2. se 4.6620 
e922 11. De 33 4, 1 2. 3 4.7313 
94 Il. 2. oe a, l Zs a. 4.7985 
.96 jl. Dene 3. 4, 1 2s 3. 4.8629 
98 jl. 2. Sr 4, I ie ce 4.9305 
1.00 [2. Ss 4, 3% 2 Se 4. 5.0000 
r d/r) ratio at nozzle throat 
-20y 1.7812 | 2.3994 : : : ; , 3.3181 
25y,, 1.7693 | 2.2978 : : ; : . 3.0826 
.33y_ 41.6798 | 2.1550 : : : : : 2.7743 





Radii of circular loci in terms of r, on which magnetic pressure equals 
hydrodynamic pressure existing at nozzle y-coordinate ratio of (f/y .); 
for various values of current and Mach Number, % = 1.30 
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TABLE II 











ontinued 

E, i eee Men l5S eee M = 20 
7o{ -21¢ |. 31c. | P_41c"s) 5 tee | Te i Teel rey |e Tee 
~40 ~u9e .046 LY exe ~ LOU ~O 0 Pe i U 
42 11.1342 | 1.7014] 2.2685]; 2.8356] 1.1357] 1.7035) 2. 2.8391 
44/1. 1275381 2. Ze 1.17051 12.7553) 2. 2.9264 
-46 ) 1. 1.8070} 2. Ss 1.2052} 1.8078) 2. 3.0130 
.48 |] 1. 1.8574] 2. 3c 1.2395] 1.8592] 2. 3.0987 
moO | 1. 1.9062! 2. on 1.2720} 1.9080] 2. 3.1801 
m2 1. 1.9558] 2. 3% 1.3050] 1.9575] 2. 3.2624 
moo | L 2.0048 | 2. 3; 1.3377] 2.0066 | 2. 3.3443 
Foo |) 1. 2.0524 | 2. 32 1.3693} 2.0540} 2. 3.4233 
58 |-1. 2.1004} 2. cP 1.4014] 2.1022) 2. 3.5036 
60] 1. 2.1482 | 2. oe 1.4333} 2.1499; 2. 3.5832 
Oz. | L.. 2.1949 | 2. 34 1.4642} 2.1964); 2. 3.6606 
64} 1. Deceit ellie. 3. 1.4951} 2.2427) 2. 3.7378 
.66 | 1. 2.2866 | 3. i 1.5252| 2.2878) 3. 3.8130 
68} 1. 223392003. 36 1.5548} 2.3323) 3. 3.8871 
MO 11. 2.3762] 3. 5S. 1.5849} 2.3774} 3. 3.9623 
my 2 | ke 2.4197 | 3. 4. 1.6139} 2.4209] 3. 4.0349 
7411. 2.4622 | 3. 4. 1.6431} 2.4647} 3. 4.1078 
760) 1. Zeo02 7 loo. 4, 1.6717) 2.5075 | 3. 4.1792 
7: an ae 2.5499 | 3. 4. 1.7004] 2.5506; 3. 4.2510 
80/1. 2.5918 } 3. 4, 1.7283 | 2.5925 | 3. 4.3208 
82,1. 2.6348 } 3. 4. L.75722 230357 \ 3. 4.3929 
84 } 1. 2360752 3% 4. 1.7840} 2.6759 | 3. 4.4599 
.86 ), 1. 2.7174 | 3. 4. 1.8119] 2.7179] 3. 4.5298 
B88 i 1. 2.7584 | 3. 4. ly. 8392).-2-7589' | 3. 4.5981 
.90] 1. 2.7982 | 3. 4. 1.8657] 2.7985] 3. 4.6642 
moe i 1. 26S 901 oe 4. 1.8934| 2.8401] 3. 4.7335 
94), 1. 2.8799 | 3. 4. 1.9201] 2.8802} 3. 4.8003 
.96 | 1. 2.918 30193.. 4. 1.9457] 2.9186 | 3. 4.8643 
.98 | 1. 2.9583 | 3. 4. 1.9722] 2.9583] 3. 4.9305 

1.00 |} 2. 3.000 4. owe 2.0000 | 3.0000); 4. 5.0000 


ratio at nozzle throat 





3.3244 
3.0890 


2.7800 


Radii of circular loci in terms of r, on which magnetic pressure equals 
hydrodynamic pressure existing at nozzle y-coordinate ratio of (f/y ), 
for various values of current and Mach Number, # = 1.30 
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TABLE II 





Continued 
M = 25 
Yo 2Ic |  3Ic | ote 
.40 | 1.100 ie Paez? 
5 a ae Ae. 2.8409 
44 11. ls 2.9303 
ma6. | 1; l. 3.0146 
4811. l. 3.1004 
250>( 15 Li. 3.1817 
eae i i. 3.2639 
5441. Zi 3.3460 
SO ulel. ee 3.4247 
ie, tal ia 3.5048 
60) 1. Ze 3.5840 
sOZe i 1. 2. 3.6617 
.64 1 1. ae 3.7389 
.66 | 1. Tage 3.8142 
68 11. ae 3.8833 
JO 41. Ze 3.9635 
Sy ae : 4.0358 
74 | 1. 2: 4.1088 
76 V1. 2? 4.1799 
781. Ze Sas) | 
80 { 1. Z SoZ) 5 
SZ. 2 4.3937 
.84 11. 2 4.4603 
200 4 1 2 4.5306 
.88 11. 4.5985 
90 I 1. Pa 4.6646 
EO2 le 2 4.7340 
94 11. 2 4.8003 
96 1. 2 4.8643 
98 J 1. 2 4.9305 
1.00 3 5.0000 


ratio at nozzle throat 
2.4150 | 2.9181 





Radii of circular loci in terms of r, 

on which magnetic pressure equals 

hydrodynamic pressure existing at 

nozzle y-coordinate ratio of (f/y_), 

for various values of current and 
Mach Number, 4 = 1.30 
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TABLE III 


é M_= 10 
Pell We | Ble | lc | SWere, 21¢ mio nema em meno 


.6727 
7633 
.8527 
9401 
.0276 
-1141 
1977 
2804 
»3652 
-4440 
.5278 
6075 
.6862 
7642 
. 8426 
.9188 
9946 
.0710 
~1442 
2205 
2915 
.3657 
4366 
. 5086 
.5813 
6494 
7210 
7916 
.8581 
9276 
0000 


NO OO el 
@ @  @  @ @ j$@ i j$@ @ @ @  @  @ 


Ze 
2. 
Za. 
2. 
25 
oe 
2. 
2 
c 
Ise 
Zz 
Ze 
Ue 
Ze 
3. 
She 
3. 
oye 
3. 
3. 
3% 
oe 
ae 
oN 
ae 
3. 
ae 
3. 
32 
3 

4 


& WD WH WD WH WD & WH WD WD WD WD W& WD W WD W W PD PR ND A NP AD PHP PNP NP AP LPL PO fo 
WMWerrrrre rPPe Pk FW W W WH WH WH WH WD W WH WH WH W PD ND NP NP 


MWe EEE SS WWW WW WD WW WW WWD DHE DP PO 
eo 
WNNNNNNNN NN NNN NN NNN NNN ND BR RRP eee 


WN NM NH NH NH KN LK NH NP NH NH KN KN PL ho 





r d/r) ratio at nozzle throat 
-20y Med 7521 2.35657 2.8713 | 3.2629 |] 1.7803 | 2.3977 | 2.8914 | 3.2894 


Radii of circular loci in terms of r, on which magnetic pressure equals 
hydrodynamic pressure existing at nozzle y-coordinate ratio of (f/y)s 
for various values of current and Mach Number, & = 1.35 
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TABLE ILI 
Continued 
M = 20 


a : [aie [ae | te [ae ie Tare ar oe) 
40 


yr. POLO) "22 2 x6 Ne 1.6126 | 2.1 2.6877 
| 1.6647 j 2. a 24 i. ee 1.6669 | 2. 2.7782 
ie 1.7182 | 2.2910] 2. 1.1469 | 1.7204 | 2. 2.8673 
| 1.7706 | 2.3609 | 2. 1.1818] 1.7727 | 72, 2.9546 
| 1.8223 | 2.4297 | 3. 1.2162; 1.8243 | 2. 3.0404 
he 1.8738 |} 2.4984 | 3. 122503 | .1..8/555) 2. 3.1258 
Ibe 1.9246 | 2.5661 | 3. 1.2841 | 1.9262 |} 2. 3.2103 
Ls 1.9734 } 2.6312 | 3. P.3167 | 1.9751 2. See lg 
ite 2.0242 | 2.6989 | 3. 1.3505 | 2.0258 | 2. 3.3763 
l. 2.0713 | 2.7617] 3. 1.3819 | 2.0728 | 2. 3.4547 
Ve ZAZIS |-2.826/ 1 5. 1.4152) 2.1228 2. 3.5381 
|e 2.1690 | 2.8921 | 3. 1.4470 | 2.1705 | 2. 3.6174 
i 2.2160 | 2.9546 | 3. 1.4782 | 2.2173 | 2. 3.6955 
| 2.2624 | 3.0166 | 3. 1.5088 | 2.2631} 3. 3.7719 
he 2.3091 | 3.0788 | 3. 1.5402 | 2.3104 | 3. 3.8506 
Le 2.30400 3 .1395:)03% 145707 52235561 | 3. 3.9268 
Le: 2.3998 | 3.1997 | 3. 1.6006 | 2.4010} 3. 4.0016 
Ve 2.4458 | 3.2610 | 4. 1.6312 | 2.4468 | 3. 4.0780 
Ll 2.4892 | 3.3189 | 4. 1.6602 | 2.4902 | 3. 4.1504 
Le 2.5349 | 3.3798 | 4. 1 .69035)° 275355 }3. 4.2258 
l. 2.5771 | 3.4361 | 4. 1.7187 | 2.5780 |} 3. 4.2966 
l. 2.6215 | 3.4953 | 4. 1.7483 | 2.6224 | 3. 4.3706 
le: 2.6636 | 3.5515 |} 4. 1.7764 | 2.6645 } 3. 4.4409 
te 2.7066 | 3.6088 | 4. 1.8049 | 2.7073 | 3. 4.5122 
1. 2.7500 | 3.6667 | 4. 1.8337 | 2.7505 | 3. 4.5842 
i 2.7910 | 3.7213 | 4. 1.8608 |; 2.7912 | 3. 4.6520 
i 2.8337 | 3.7782 | 4. 1.8893 |; 2.8339 | 3. 4.7232 
1. 2.8755 | 3.8340 | 4. 1.9172 | 2.8758 | 3. 4.7929 
We 2.9155 | 3.8873 | 4. 1.9436 | 2.9154 ; 3. 4.8591 
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ratio at nozzle throat 


1.7811 | 2.3999 | 2.8951 | 3.2920} 1.7816 | 2.4007 | 2.8962 | 3.2956 


Radii of circular loci in terms of r, on which magnetic pressure equals 
hydrodynamic pressure existing at nozzle y-coordinate ratio of (f/y, ); 
for various values of current and Mach Number, ®& = 1.35 
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TABLE III 

(Continued) 

M = 25 
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r d/r) ratio at nozzle throat 
-20y , 1.7817 | 2.4010 | 2.8969 |} 3.2964 


Radii of circular loci in terms of r, 

on which magnetic pressure equals 

hydrodynamic pressure existing at 

nozzle y-coordinate ratio of (f/y_), 

for various values of current and 
Mach Number, x = 1,35 
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Radii of circular loci in terms of r, on which magnetic pressure equals 
hydrodynamic pressure existing at nozzle y-coordinate ratio of (f/y, v7 
for various values of current and Mach Number, & = 1.40 
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TABLE IV 
xCone tenes 
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1.7758 | 2.3875 } 2.8751] 3.2689 


Radii of circular loci in terms of r, on which magnetic pressure equals 
hydrodynamic pressure existing at nozzle y-coordinate ratio of (f/y ); 
for various values of current and Mach Number, % = 1.40 
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TABLE IV 


(Continued) 
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r (d/r) ratio at nozzle throat 


. 20y . 1.7760 | 2.3879 | 2.8757 | 3.2696 


Radii of circular loci in terms of r, 

on which magnetic pressure equals 

hydrodynamic pressure existing at 

nozzle y-coordinate ratio of (f/y_), 

for various values of current and 
Mach Number, % = 1.40 
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Ce aes M = 10 
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d/r) ratio at nozzle throat 
' 
1.76371 2.3582 | 2.8289 | 3.20908 1.7690 | 2.3711 | 2.8499 


Radii of circular loci in terms of r, on which magnetic pressure equals 
hydrodynamic pressure existing at nozzle y-coordinate ratio of (f/y.), 
for various values of current and Mach Number, @ = 1.45 
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TABLE V 
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nozzle throat 


}1.7700 | 2.3735 2.8537 3.2613] 1.7703 | 2.3743 





3.2431 


Radii of circular loci in terms of r, on which magnetic pressure equals 
hydrodynamic pressure existing at nozzle y-coordinate ratio of (f/y)> 
for various values of current and Mach Number, @ = 1.45 
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TABLE V 


(Continued ) 
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Radii of circular loci in terms of r, 

on which magnetic pressure equals 

hydrodynamic pressure existing at 

nozzle y-coordinate ratio of (f/y_), 

for various values of current and 
Mach Number, @ = 1.45 
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TABLE VI 
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ratio at nozzle throat 
1.7578 12.3452 | 2.8089 | 3.1813 ]2.7635 2.3582 | 2.8298 | 3.2104 


Radii of circular loci in terms of r, on which magnetic pressure equals 
hydrodynamic pressure existing at nozzle y-coordinate ratio of (f/y)s 
for various values of current and Mach Number, fY = 1.50 
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TABLE VI 
Continued 
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nozzle throat 
1.7648 | 2.3613 






1.7645 | 2.3605 





2.0398 goe 2197 2.8347 


Radii of circular loci in terms of r, on which magnetic pressure equals 
hydrodynamic pressure existing at nozzle y-coordinate ratio of (f/y), 
for various values of current and Mach Number, 2 = 1.50 
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d/r) ratio at nozzle throat 
.20y 11.7650 | 2.3617 | 2.8353 13.2178 
O 


Radii of circular loci in terms of r, 

on which magnetic pressure equals 

hydrodynamic pressure existing at 

nozzle y-coordinate ratio of (f/y _), 

for various values of current and 
Mach Number, @ = 1.50 
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TABLE VII 


CROSS-INTERPOLATION TABLE 
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Perpendicular distance in terms 
of r from current-carrying 
conductor to lines parallel to 
channel wall within the channel 
at ordinate f/y. 
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